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ABSTRACT 


The problem of determining the limiting performance characteris- 
tics of mechanical systems subject to random input is studied. A 
review is presented of the classical work in the optimal design of stochas- 
tic systems. Some recent results of stochastic optimal control theory 
are employed. The solution to the limiting performance problem is 
formulated in both the frequency and time domains. Both formulations 
require substantial, burdensome computations when applied to large scale 


systems . 



1* INTRODUCTION 

Random disturbances appear as complicated time-varying functions 
that may exhibit wide irregular variations in amplitude and frequency. 
Both the input disturbance and the system response must be given sta- 
tistical characterization and. as expected, this complicates the optimum 
design problem. No all encompassing methodologies are available for opti 
mlztng realistic systems under general random environments. Since 
related literature from control theory on the optimization of stochastic 
processes is quite extensive, we will take advantage of these develop- 
ments and use these techniques, both in the frequency domain and in the 
time domain, as a basis for this study of limiting performance of dynamic 
systems subject to certain classes of random inputs. The limiting per- 
formance concept for transient systems is defined in Ref. 1. Briefly, 

it is the theoretically optimum performance of a system, regardless 
of configuration. 

In this study the input disturbances are characterized by spec- 
tral density matrices or correlation matrices. Although for much 
of this report the performance index is quadratic in both responses and 
control forces, other choices are possible. Most of this study is con- 
fined to linear systems. 

This Study begins with a review of the single-input , single-output 
Kiener-Hopf type approach which has been examined by others (Ref. 2,3). 

The governing Wiener-Hopf equation and its solution are summarized. An 
example is presented. Next, it is shown how the limiting performance of 
a multi-degree of freedom (MDF) system under random input can be deter- 
mined by the extension of Wiener-Hopf techniques to the multiple-input, 
multiple-output case reported by Weston, et. al., (Ref. 4). 
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Parallel to the above frequency domain appreachj the problem in 
the time domain will also be formulated by applying results from stoch- 
astic optimal control theory (Ref. 5). This approach has also been 
pointed out recently by Karnopp (Ref. 6). a single degree of freedom 
(SDF) problem which was previously solved by frequency domain techni- 
ques will be solved by this approach. 

Also, included in this report will be dynamic programming ap- 
proaches. This method is more general In the sense that it can treat 
performance indexes other than a mean-square type. A formulation with 
quadratic criteria is summarized and a non-quadratic possibility pro- 
posed. 

1.1 Statement of Problem 

Consider a MDF dynamical system described by 

^ + jCx + ^ + Vu « Ff (1) 

where 

M - n X n mass matrix 
£ - n X n damping matrix 
K - n X n stiffness matrix 

V - n X nu coefficient matrix of control forces 

JF - n X nf coefficient matrix of disturbances 

X - n X 1 state vector 
f - nf X 1 disturbance vector 

u - nu X 1 generic or control force vector 

n - number of degree-of-freedom of the system 
nu - number of controllers of the system 

number of disturbances (forcing functions) of the system 


nf - 
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In Eq, (1) ? is a random vector of known statistics* In order to make 
a limiting performance study, portions of the system have been replaced 
by u a vector of generic or control forces These forces can represent 
any sub— system. The remainder of the system must be linear as must be 
the system kinematics. 

Since the forcing function is random in nature, the limiting per- 
formance problem must be defined in some statistical sense. The most 
common optimization criterion is to minimize some expected mean square 
response quantity while imposing constraints on other expected mean 
square responses. This criterion, though not as direct as the criterion 
used in the transient case is still meaningful. It is selected because 
in the analysis of stochastic processes, it is difficult to find responses 
other than the mean square response of the system. Another reason for 
using mean square criteria is that when the performance index and con- 
straints are combined in a penalty function type objective function, 
the resulting optimum system will be linear* This is a well known 
result of optimal control theory (Ref. 5). Other criteria have been 
suggested by several authors. Aoki (Ref. 7) treated a simple problem 
based on a performance index of a maximum expected deviation over the 
history of a dynamic process. Trikha and Karnopp (Ref. 8) studied a 
single degree of freedom system with a criterion based on the values of 
displacement and acceleration, whose probability of being exceeded is less 
than a prescribed value. It does not appear as though either of these 
criteria can be applied to large. systems. 



2. REVIEW OF CLASSICAL WO^ 

Limiting performance characteristics of SDF systems have been 
reported in the literature for several inputs and are based on either 
expected mean square values or the probability of exceeding selected 
response levels. These results are summarized in Ref. 1. To make this 
report more self-contained, the limiting performance problem based on 
expected mean square responses will be treated briefly. 

2.1 The Governing Integral Equation 

Consider the configuration of Fig. 1, For a given input spectral 
density the optimum compensation W^(t) (or W^(a) in the frequency 

domain) is to be found such that E[y^(t)] is minimized subject to the 
constraint ]^o, where a is the maximum allowable mean square value 

and E[’] stands for mathematical expectation. This constrained 
optimization problem will be converted to an unconstrained problem by 
introducing a Lagrangian multiplier p. The problem then becomes one of 
minimizing 

F = E[y^] + pE|‘q^} (2) 

It is shown in Ref, 9 that the optimum compensation W (t) that 

cm 

minimizes (2) is governed by the integral equation 
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v(t) = random input disturbance 

q(t) = actual output 

i(t) = ideal output 

y^(t) = error 

Fig. 1 Block diagram for deriving Wiener-Hopf 
integral equation. 
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where 


r<ti) - (tj^+t2)dt. 


° t/” r8f(t2)gj(t^)dt2dt^ 


«CQ ^ ^OO 




= cross correlation function of input to ideal output 


(ft = autocorrelation function of input 

w 


2.2 Solution of the Governing Integral Equation 

Equation (3) is of the Wiener-Hopf type for which spectral fac- 
torization can be applied to obtain an explicit solution formula in the 
frequency domain (Ref. 9). This is 


cm 


r(s) 


A_LsU+ 


+ a 


A^(s) 


(4) 


where a Is an arbitrary constraint 


r(s) = 2irG^(-s) 


A(s) = 27 t [G,(s)G,(-s) + pG„(s)G (-s)] $,„X8) 


A (s) ~ any factor of A<-a) which includes all the poles and zeros 
in the LHP (left-half-plane) 


A (s) = A(s)/A (s) 


IIsl 

_A“(s)J 




- component of which has all its 


poles in the LHP such that 


r(s) 

A’'(s) 


r(s) 


A- (s) 

L- _l + 


has all its poles 


in the RHP , 
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2-3 Application to SDF System 

The results of previous sections can be applied to find the limiting 
performance characteristics of a single degree system subject to random 
inputs. Consider the SDF system shown in Fig. 2. The problem is to find 
u such that the performance index E[u^J + pE[x^] is minimized.- By 
the present approach, we seek to find the optimum transfer function 

\j(a) = 

Draw the block diagram of Fig. 3. Equation (4) will provide W(s) for 
given input spectral density. 

For an input with spectral density 


The solution is 

W(s) - e^/(s^ + /I s6 + B^) 

where 



2.4 Discussion 

The classical Wiener-Hopf technique leads to tractible analytical 
solutions to single- input single-output systems. Since most systems do 
not fit in this class the usefulness of this technique is quite limited. 



y(t) (input motion) 


Fig. 2 SDF system subject to random input. 
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G^(s) = ms^ Gf(s) = 1 H(s) = 1 

W(s) * ? 

Fig. 3 Block diagram for the configuration of Fig. 2. 
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3. EXTENSION OF WIENER-HOPF TYPE APPROACH TO THE LIMITING PERFORMANCE 
OF MDF SYSTEMS 

The analytical design techniques of Newton et;al,(Ref, 9), have 
been applied to find the limiting performance characteristics of simple 
systems. These techniques are extended by Weston and Bongiorno (Ref. 4). 
It is the purpose of this section to summarize the principal results of 
this extended theory. Also, an example will be given to Illustrate the 
application of this technique. 

3,1 Formulation 

Consider the system configuration shown in Fig. 4. All elements 
of the system are assumed; to be linear and time Invariant. The input 
vectors (signal, noise, and disturbance) are generated by a stationary 
stochastic process whose power spectral density is known. The purpose 
is to find optimum compensation to minimize a performance index con- 
sisting of a weighted sum of the output mean-square error plus a weighted 
sum of the mean square value of a set of saturation signals. The results 
are summarized in Theorem 1. 

The following nomenclature is used. - r^(t), n(t) , d(t) are n x 1 
vectors of input, noise and disturbance (to the measurements) respec- 
tively, whose elements are the realization of an independent, stationary 
random process* The statistics of the random processes r^^, n, d are 
assumed to be known and adequately described by rational power spectral- 
density matrices „ (s) , ^j^(s), ^o(s), respectively. Transfer 
function matrices ^(s), ^(s) and ll(s), represent asymptotically stable 
systems. The dimensions of these matrices are; 

^-n X m (m^n), ^ -m x n, H,-h ^ -q x m 

q is the dimension of A(s) , the satuation signals. 



R^(s)- 

input 

signal 


/ 
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N(s) 


input 



disturbance 
D(s) 


) — 

G (s) 


G (s) 
P 

i 

f 

c 





compensator 


feedback 


plant 


actual system 
Ideal system 






error 

Y (s) 
e 


Fig. 4 Block diagram of multivariable system. 
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The performance index to be minimized is 




( 5 ) 


where ^ and R are real, constant, S 3 nnmetrlc, positive, definite matrices. 

£isnxn,Ris qxq. The superscript T designates the transpose of a matrix. 

The problem was solved by a calculus of yariatipn .technique with the results 
summarized in Theorem 1. 

Theorem 1 . (Ref. 4) 

The physically realizable W - W' that minimizes 

c cO 

J = 

associated with the system in Fig, 4, satisfies 


I«c0 “ [#]"’■ 


( 6 ) 


where 




1*1 “ Gp*S!Gp + 2s*^% 


a) 


A = \l5Gp* 


The optimal compensation is then 


G = [I-W _ HG W ^ 

~c •~c0 — •’• p -*c0 


( 8 ) 


The following notatlonal convention is used in Theorem 1; 


(1) For real matrix X, 2^(s) = 2T(-s) 


(2) [ZCs)]^ = part of function Z(s) which is analytic in R^s^O 
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3.2 Application 

For any nmltl^inputi multi-output system, the application of 
Theorem 1 leads to the optimum transfer matrix which minimizes the per- 
formance index of (5), provided a block diagram similar to Fig. 4 can 
be constructed for this system. We will illustrate this using a two 
degree of freedom system. This same problem was solved by Bender 
(Ref. 3) using the classical Wiener-Hopf approach. 


Consider the system shown in Fig. 5. The input to this system 
has the following spectral density 


$ (s) = 

XX 2 

s 


This system is governed by the differential equations 
My = u 

mi = -u + k(y-z) 


( 9 ) 


- The problem. Is to find the optimum isolator characteristics such 

1 ' ' V' ' 

that .the relative displacement (z-y) la minimized while the acceleration 
z is bounded. Since the input is random, a meaningful performance index 
will be of the mean square type 


J “ E[py^^+ (z-y) ^3 


The problem will be solved in the frequency domain- 
Define , 

W(s) * Y(s)/XCs) 

5(s) Z(s) - Y(b) 


( 10 ) 


Cll) 

( 12 ) 
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Transforming (9) into the frequency domain it is found that 


Ms Y(s) = U(s) 


ms Z(s) = -U(s) + K(X(s) - Z(s)) 


Solving (13) using definitions (11), 


(13) 


Z(s) 

K+s m 


(14) 


Y(s) - X(s) 


(15) 


6(s) « Z(s) - Y(s) 


'm 


[. (1+1/r) (s/w )*^ + 1] 
n •* 

s^(l+s/tp )^ 
n 


W(s) 


X(s) 


where = k/m, r = m/n 


From the block diagram of Fig, 6 


^(s) * [H 2 (s) - H^(s) K(s)] x(s) 


Compare the above two equations, 


H,(s) 2 


1 + (1 + 1/r) (s/o. ) 


n 


s‘(l + (sMjj)^) 


(16) 
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vibration 



Fig. 6 Block diagram for Fig. 5. 
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H,(s) = i J (17) 

1 + (s/u^)‘‘ 

Compare the block diagrams of Fig. 6 and Fig, 4. For the present 
case, the various matrices used in Theorem 1 are 

H = 0, G = G = H- , G = W, G =1 
i 2 p 1 c s 

All quantities are scalars for this system. For the performance index, 

Q - 1, R = p. From Theorem 1, we have for this case 

"cO = (18) 

where 

= -kis (19) 

V-Vi+p 

A » H2Hj^j^(-A/s^) (20) 

From (19) f! = Tkis (21) 

Substitution of (19) , (20) , (21) , into (18) gives the same equation as 
found in Ref, 3. 

3.3 Discussion 

It would appear that the results given in Theorem 1 are quite prom- 
ising, in that multi- input, multi-output systems can be treated. However, 
there are at least two difficulties involved. First of all, to use 
Theorem 1, a block diagram similar to Fig, 4 must be established for the 
system. Unfortunately, this would be a difficult procedure to automate. 
Secondly, computational problems are encountered when this is extended 
to large systems. As indicated in Theorem 1, the solution of the optimum 
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transfer function matrix requires matrix spectral factorization. For 
practical problems, this must be d on e computationally. An. .available pro- 
gram using Tuel’s algorithm (Ref. 10) can be used to factor polynomial 
matrices of dimension up to 6x6 with order not exceeding 20 for each 
polynomial element. This program requires 10,000 words of storage. 
Apparently, this is a very limited computational capability, 

4. APPLICATION OF STOCHASTIC OPTIMAL CONTROL THEORY TO THE LIMITING PER- 
FORMANCE CHARACTERISTICS OF MDF MECHANICAL SYSTEM 

The formulations described in the previous two sections are fre- 
quency domain approaches. In this section a time domain approach will 
be set forth, with the equations Qf motion in a form compatible with those 
used in the limiting performance study of transient systems (Ref. 11) 

4.1 General Formulation 

Consider a mechanical system described by Eqs. (1) with a 
gausslan white noise representation for f. Consider a linear realiza- 
ble system for u. Define 

Rj(t) = E[I (t)f^(t+T)] 

as the correlation matrix of the input disturbance vector T. For this 
problem 

^(t) « ^6(t) (22) 

where 6 is the usual Kronecker delta function. 

Since the input disturbances are random in nature, the limiting 
performance characteristics will be defined in the expected mean square 
sense. That is, we will find u that minimizes the expected mean square 



of some response (objective function) while the expected mean square values 
of some othe r re sponse variables (constraints) are bounded. Trade-off 
relations between the objective function and any one of the constraints 


can then be obtained. 


To solve this problem, advantage will be taken of stochastic opti- 
mal control theory. To do so we transform our problem to the format 
of optimal control theory, A similar technique has been recently 
reported by Karnopp (Ref. 6). 

Equations (1) can be written in first order form as 


s “ As + Bu + Gf 


where s - jx 

k 


A * 


-1 -1 
-M -M C 


( 23 ) 


B 


0 



The limiting performance problem is to find u(t) such that the performance 


index 
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is minimized. In (24) E(<*] stands for expectation. ^ and R are 
positive semidefinite weighting matrices* By adjusting the components 
of the ^ and R matrices properly, the desired limiting performance 
problem can be formulated.; 

In stochastic control theory, an observation (or measurement) 
model is required. Since we are considering the limiting or theoreti— 
cally optimum performance, it will be appropriate for us to assume 
noise free perfect measurements. For such problems, optimal control 
theory (Ref. 5) leads to the following results. For the system des- 
cribed by (23) Subject to white noise input, the optimal control law 
u* that minimizes the performance index (11) is given by 

U* = (25) 

where Y is a constant matrix which is the solution of the following 
algebraic matrix Riccati equation 

- YBR^^B^Y + ^ = 0 (26) 

^^ving the optimal control law, the next step is to calculate the 
value of the performance index and then establish the desired trade-off 
relations. To evaluate the correlation matrices such as E[xx'*’], a 
method reported by Yang and Iwan (Ref. 12) will be applied. To do so, 
substitute the control law into Eq, (23) and rewrite the resulting equa- 
1^ second order form. Then, by application of the result given in 
Ref. 11, the following set of linear equations for the computation of 

X --T 

E[xx ],’E[xx ], and E[xx ] are obtained. 
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T 


I 

+ c^ E[kx^]M^ + ^ [xx^jK^ + K = fR.F^ 


—T. 


(27) 


where 


£q 


C+JY22 


■ K - K + TY_ 

0 — 21 

T * VR~V^» 


Y Y 

“11 -12 


Y Y 

-21 -22 


and 


E[uu^] * Z^ELx^’^lZ^ + 


(28) 


-1, 


where [^2? ^ 


^•2 Example: Limiting Perfortaance of an SDF Vibration Isolater 

Consider the SDF system shown in Fig, 7. Suppose the disturbance 
f is gaussian white noise. The problem is to find the trade-off rela- 
tion between mean square values of k ( or u/m) and x. This problem has 
been solved by Fujiwaira et. al. (Ref. 14). We will solve the problem 


in the time domain. 
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The governing equation can be written 


as 


nuc = u + f 


(29) 


or in first order form 


s =-As + Bu + Gf 


(30) 


where 


r ^ 1 

8 = { . } 


- RjS(t) 


A » 


0 1 
0 0 


B = 


0 

1/ra 


G = 


- Li% 


The performance index to be minimized is 


J = 


E[/^ + s"^Qs) dt] 

0 


(31) 


where 


R = il] 


0 

0 0 


The optimum control law is 


-1 T - 
u* » -R B Ys 


(32) 


where Y is the solution of 


A^Y + YA -* YBR”^B^Y + Q = 0 


(33) 
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The solution of (33) is 



A /2mA 

mA 

I == 

mA 

m/2mA 


Substitute (34) into (32) to find 

u* = '^s^(t) - s^(t) 

or 

u*(t) == -Xx(t) - /2mA x(t) (35) 

Since E[xx J = 0 for an SDF system, (27) becomes 
ME[x^] - K E:[x^] = 0 

ME[i^]Cp + CgE[x^]M + 0 - (36) 

where 

Kq “ a, Cq = 

Thus. E[x^] - R^/2A/^ and 

E[u^] = -.(3 r^/2.^) /?T 
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These results are the same as obtained by Fujlwara. By eliminating ^ 

2 2 

between E[u ] and E(x ] , we find the tradeoff relationship 

(^) * 7 

where 

b = R,/2vS 
f 

This is plotted in Fig, 8. 

4.3 Discussion 

The underlying assumption for the formulation of this section is 
that 1 must be white noise. This is not so severe as it may appear, even 
though white noise does not actually exist. White noise does provide 
a good approximation for broad band random noises. The major drawback 
in this formulation is that its implementation encounters computational 
burdens for large systems. 

Since the most efficient procedure of solving matrix Riccati equa- 
tions (26) is by matrix spectral factorization (Ref. 4), this approach 
involves the same complications as described in Section 3. However, set- 
ting up this formulation for a particular problem is much easier then 
using the method of Section 3. 

5. DYNAMIC PROGRAMMING APPROACH 

In this section the application of dynamic programming techniques 
to the limiting performance problem of dynamic systems subject to random 
inputs will be discussed. This technique is more powerful in the sense it 
can treat problems with performance indices other than the usual expec- 
ted mean square type. The essence of this technique will be discussed 
first and then applied to various limiting performance problems of interest. 
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5.1 Definition of Dynamic Progracnalng 

Dynamic Programming is the name given by its Inventor, Richard Bellman 
(Ref. 15), to a computationally motivated procedure for solving optimi-" 
Eation problems through a sequence of smaller problems. It is ideal for 
multi-stage optimization problems since it provides an efficient algo- 
rithm for analytical and computational procedures by allowing only those 
continuations of multi-stage processes that constitute optimal continu- 
ations. This can be formalized as the Principle of Opf -i t-y , This 
principle can be applied to derive the functional relationships of dy- 
namic programming. Many researchers have applied the dynamic programming 
technique to solve stochastic optimal control problems. A rather com- 
prehensive bibliography in this regard can be found in Ref, 16. 


To illustrate the principle of optimality, consider the following 
simple example: Given: 


®i+l T(Si, u^) 


(37) 


with s^ = c where s =* state vector, u = control vector 


find • ‘ » •» ^-1 minimize 


N 


't' = I Hs.) 
i=l 


(38) 


Define 


fj^(c) = Min Min . 

^0 ^1 


. Mln[(J)(s^)+ 


. + 


(39) 


u 


N-1 
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Then, by definition, 

f (c) - Min 4>(T(c, u )) (40) 

“0 

for N12 

f^(c) - Min [4 >(t(c,Uq) +%i(T(c,Jq)] (41) 

The principal of optimality, when applied to the above example, 
specifies that when only one more decision (u’s) stage remains, the 
decision must be made in such a way as to minimize <t> (s^) . When more 
than one decision stage remains, make the immediate decision in such 
a way as to minimize the overall sum consisting of the terms of the 
immediate contribution and the term fj^_^{X(c,UQ))^ representing sums 

obtainable from the resulting state T(c,Uq) by an optimal continuation* 
Solution of the functional relations (41) gives the optimal control 
sequences . 

In the following subsections, several limiting performance problems 
will be formulated with different optimization criteria. These criteria 
include: mean quadratic performance index, probability of maximum deviation* 

5.2\ Quadratic. Performance Index 

' . ; ■ . II tF 

Consider a mechanical system described by (23) in which T, the 
disturbance, is an independent random variable. The discrete solution 
of (23) can be written as 

s(k+l) *= $s(k) + Gu(k) + h(k) 


(42) 



where 


^ ■ er- T = sampling period 

G = $(t-T)B dx 

kT — ~ 

h(k) = J^(t-T)D?(x)dT 

kT 

The performance Index of Interest is the mean quadratic form 
N -T 

^ (k)£s(k)'+y^(k-l)Ru(k-l)] ( 43 ) 

k -1 

where £, R are positive definite symmetric matrices. By adjusting the 
values in the elements of ^ and R matrices, relative weight can he 
Imposed on state as well as control. Following a dynamic programming 
approach as used by Tou (Ref. 17), the following results can be derived. 

The optimal control law that minimizes the performance index (43) 
can be obtained from the algorithm 
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where 

K(K-J) = -(Lgg(N-(j+l)+R)"^^^(N-(j+l)) 

I (N-(j+D) = I^^(N-Cj+1)) +L^g(N (J+1 ))B(n-J) 

l^^(N-(j+l)) ■ l'^S(N-(j+l))4.(J) (45) 

icG^N-O+D) = G’'^S(N-(j+l))G 

ig^W-O+l)) = G\(N-(j+l))* 

L)G(N-(j+D) = $^s(N-(j+l))G 

S^(N-(j+l)) =2.+ P(N-(j+l)) 

with 

P(0) = 0 

Thus, from j « N - 1, with _P(0) = 0 we can generate K(l) , P (1) , K(2) , 
£(2), . . . recursively. Note that even when the system is invariant, 
the feedback gain matrix Is still tlaa varying. The performance index 
can be evaluated from 

f^[x(0)] = x^(0)P(N)x(0) 

Such problems as the airplane taxiing during landing can be treated 
by this technique. The landing velocity provides an initial condition 
and the runway is the random Input. Also, taxiing stops after a finite 
time. This approach can be used to evaluate the limiting performance 
for such systems. 
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MJ-PMilzlnR the Prob ability of a Maximum Devlatlnn 
Bellmai, (Ref. 18) pointed out that the problem of minimizing the 
probability that the maximum response will exceed a preassigned value z, 
can be solved by dynamic programming. For this case 

£^(c) “ Prob{max(|^| . . .. |^|)<i) 

(46) 

where |x|£z 

If we interprete |x|<z. in the component by component sense, then 
this problem, when solved for a range of z, will result in an interesting 
trade-off diagram (Fig. 9). For example, for z = A in Fig. 9, the 
probability that max |xj will not exceed A will be . The probability 

ofe exceeding A is (IJ?^). This could be correlated to some other prop- 
erties of the system, such as reliability. 

5.4 Discussion 

Dynamic programming provides a very general limiting performance 
problem formulatlnn. However, due to its Inherent "curse of dimensionality"^ 
and the concomitant computational burden, the application of this techni- 
que to large scale problems is limited. 

6. CONCLUSIONS 

The limiting performance characteristics for dynamic, systems sub- 
ject to random Inputs has been solved for a class of problem by various 
control theory optimization techniques. In most cases, these problems 
are restricted to linear systems subject to Inputs with known power spec- 
tral densities, either white noise or representable by rational functions. 
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In the case of white noise, stochastic optimal control theory in the 
time— domain can be applied to find the optimal law leading to limiting 
performance. This approach has been formalized for vibration isola- 
tion systems and only requires matrix manipulation and solution of a 
set of algebraic Riccati equations. 

The analytic design technique extended by Weston is more general 
in the sense that it can solve problems subject to colored noise inputs. 
The drawback in this approach is that the system configuration must be 
put in a special block-diagram form. However, for the computation of an 
optimum transfer function, existing spectral factorization programs can 
be applied. 

Dynamic programming, which appears to be the most versatile formula- 
tion, is cf limited value due to its computational inefficiency. 
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